Theory of Tunneling in the Exciton Condensate of Bilayer Quantum Hall Systems 
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We develop a theory of interlayer tunneling in the exciton condensate of bilayer quantum Hall 
systems, which predicts strongly enhanced, but finite, tunneling conductance peaks near zero bias 
even at zero temperature. It is emphasized that, though this strongly enhanced tunneling originates 
from spontaneous interlayer phase coherence, it is fundamentally not the Josephson effect. Because 
of strong interlayer correlation, the bilayer system behaves as a single system so that conventional 
tunneling theories treating two layers as independent systems are not applicable. Based on our 
theory, we compute the height of conductance peak as a function of interlayer distance, which is in 
good agreement with experiment. 

PACS numbers: 73.43.-f, 73.21.-b 
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When Spielman et al. [ij observed strongly enhanced 
interlayer conductance peaks near zero bias in bilayer 
quantum Hall systems at total filling factor iAp = 1, they 
not only renewed our interest in the bilayer quantum 
Hall effect Q, but also attracted intense interest from 
the general perspective of strongly correlated physics. It 
was because, in addition to its many-body origin, the 
bilayer quantum Hall effect bears a rather precise anal- 
ogy to superconductivity; the ground state of bilayer 
quantum Hall effect at interlayer distance d/lg <C 1 
(Ib = y/hc/eB) maps onto the BCS wavefunction of an 
exciton condensate of particle-hole pairs formed across 
the interlayer barrier. In fact, Bose-Einstein condensa- 
tion of excitons in semiconductors has been sought after 
for decades. In particular, there have been fascinating 
recent experiments on the possible condensation of opti- 
cally generated indirect excitons |3|, for which, however, 
there is not yet conclusive evidence. On the other hand, 
it is generally accepted that the strongly enhanced con- 
ductance peak in the quantum Hall regime is a direct 
indication of macroscopic phase coherence. 

To be concrete regarding the mapping between the su- 
perconductivity and bilayer quantum Hall effect, let us 
write the exact ground state wavefunction at d/lg = 0, 
i.e. the Halperin's (1,1,1) state (which is adiabatically 
connected to the ground states at sufficiently small, but 
finite d/ls)'- 
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where m is a momentum index in the lowest Landau 
level and the pseudospin representation is used: f (j) 
indicates the top (bottom) layer. Note that Eq.Q de- 
scribes the full wavefunction including both orbital and 
layer degree of freedom |||. Since Eq.Q has a struc- 
ture isomorphic to the BCS wavefunction, it is clear that 
the bilayer quantum Hall state should have a phase co- 
herence between states with different interlayer number 
difference in analogy with phase coherence between dif- 
ferent number eigenstates in superconductivity, which is 



the origin of the Josephson effect. Naturally, this simi- 
larity led previous authors 0,01 to predict the Josephson 
effect in bilayer quantum Hall systems. The strongly en- 
hanced conductance observed by Spielman et at, there- 
fore, seemed to be exactly the experimental verification 
needed. However, there are key properties of the con- 
ductance peak indicating that this phenomenon is not 
the conventional Josephson effect: most notably, satura- 
tion of height as well as width to finite values in the limit 
of zero temperature ||. 

This apparent discrepancy gave rise to two groups 
of thought. In one group, the enhanced conductance 
is still regarded as DC Josephson effect, but its height 
is reduced by complicated disorder-induced fluctuations 
HEIHII3- On the other hand, others [3 argued that 
there is no exact analog of Josephson effect in interlayer 
tunneling experiments because the bilayer system as a 
whole is a single superfiuid, not a set of two superfliud 
systems. While we agree with the latter viewpoint that 
the enhanced interlayer tunneling conductance is not the 
analog of Josephson effect, we show below that strong in- 
terlayer correlation requires a fundamentally new start- 
ing point different from all of above theories in order to 
construct a self-consistent theory of interlayer tunneling 
in quantum Hall regime. 

As mentioned previously, the bilayer quantum Hall sys- 
tem is a single superfiuid system. So, it is impossible 
to induce a chemical potential gradient between the two 
layers without destroying interlayer phase coherence, in 
which case the interlayer current becomes a normal cur- 
rent, not supercurrent. It is important to distinguish be- 
tween the chemical potential gradient and applied inter- 
layer bias voltage because, even when the bias voltage is 
applied, bilayer systems will immediately reach an equi- 
librium by creating charge imbalance in order to compen- 
sate the relative voltage difference and therefore there is 
no chemical potential gradient. Though this point seems 
straightforward, it has been completely overlooked by 
all previous theories which, regardless of their viewpoint 
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FIG. 1: Schematic diagram of tunneling measurement in bi- 
layer quantum Hall systems. Note that there is no chemical 
potential gradient between layers when the ground state of 
bilayer system becomes a single exciton condensate at small 
interlayer distance d, as depicted in (a). A consistent the- 
ory of interlayer tunneling, therefore, should inevitably take 
external leads into consideration. On the other hand, when 
d is sufficiently large as shown in (b), two layers behave as 
independent systems, and interlayer coherence is lost. 



regarding the analogy with Josephson effect, began by 
implicitly making a self-contradictory assumption that 
there is strong interlayer correlation due to the Coulomb 
interaction but two layers can be treated independently 
by having a finite chemical potential gradient. In fact, if 
one can induce a finite chemical potential gradient while 
maintaining interlayer phase coherence, there would be 
a very interesting experimental consequence: oscillating 
tunneling current whose frequency is proportional to the 
applied bias voltage. However, no oscillating current has 
been observed in experiments. 

Now, if there is no interlayer chemical potential gra- 
dient, there is no electromotive force within bilayer sys- 
tem and any current should be induced from outside. 
It is, therefore, necessary to take into account external 
leads, as schematically shown in Fig^ This, of course, 
makes any quantitative prediction dependent on the way 
in which bilayer systems are connected to external leads. 
However, it is still possible to make a quantitative predic- 
tion on essential aspects of coherent interlayer tunneling. 
In particular, we will compute the dependence of tunnel- 
ing conductance peak height on interlayer distance cL/Ir- 
Also, we will show that the width is finite even at zero 
temperature, and it is controlled ultimately by extremely 
small, but finite single-particle interlayer tunneling gap 
Asas- 

Let us begin our quantitative analysis by writing 
the total Hamiltonian including the Hamiltonian for 
Coulomb interaction between electrons in bilayer system 
Hq, the Hamiltonian describing the left and right lead, 



Hl and Hr respectively, and tunneling between leads 
and the bilayer system H': 



H = Hq + H' + Ht. 
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H' = yT R ^{k,m)[& R {k)c m -[ + h.c] 



+ ^T ii (p,m')[4(p)c mU + h.c.] ! (4) 

where, as before, the pseudospin representation is used, 
and Vlll is the lowest Landau level projection opera- 
tor. Tn^(k,m) is the tunneling amplitude between the 
state with momentum k in the right lead, and the state 
with to in the top layer of bilayer system. Tli(p, to') 
is similarly defined. Hr and Hl describe electrons in 
external leads as normal Fermi liquids. It is now very 
important to note that H does not have any interlayer 
tunneling term within the bilayer system. It is because 
we are interested in the spontaneous interlayer coherence 
which occurs in the limit of zero interlayer tunneling gap: 
AsAs/(e 2 /e^s) — ► 0. As will be shown later, this sponta- 
neous interlayer coherence is due to the many-body effect 
of Coulomb interaction in Hq, and it creates a non-zero 
current from one layer to the other even in the limit of 
zero interlayer tunneling gap (of course, in unbiased equi- 
librium, the net current is zero since two opposite cur- 
rents cancel each other). 

Since there is no direct process of transporting elec- 
trons from one lead through the bilayer system to the 
other lead, one has to consider second order tunneling 
processes: 



H^p — H - 



1 



Hq — Hr — Hl 



H', 



(5) 



where E g is the ground state energy of Hq + Hr + Hl- 
By adding an electron to the top layer and removing an- 
other from the bottom layer, H' T describes tunneling pro- 
cesses through the bilayer system. Now, because the bi- 
layer quantum Hall state is incompressible at sufficiently 
small d/ls, adding or removing electrons costs a finite en- 
ergy which is equal to the Coulomb self-energy of quasi- 
particles, Ac We will compute Ac as a function of 
cL/Ir later by using exact diagonalization. It is, however, 
sufficient at this stage to know that Ac is independent 
of momentum m in the lowest Landau level. So one can 
just replace Hq + Hr + Hl — E g by Ac- Remember that 
there is no energy cost in taking electrons from external 
leads because normal Fermi liquids are compressible. 

Now, we assume that the tunneling amplitudes 
TRf(fc, to) and TL[(p,m') are more or less independent 
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f = ^S m T fli (m)cJ„ t c m | 



FIG. 2: Feynman diagram of interlayer tunneling in bilayer 
quantum Hall systems. The vertex operator T contains all 
of many-body effects of an exciton condensate. Tul is the 
tunneling amplitude and Ac is the Coulomb self-energy of 
quasiparticle. 



of momenta k and p, which is a common practice in tun- 
neling theories when studying tunneling processes only 
within a narrow region of energy near Fermi surface. 
Keeping only terms of H' T relevant for transporting elec- 
trons from one lead to the other, we arrive at the follow- 
ing tunneling Hamiltonian: 



k.p 



[4(fc) Ci (p)ft + c[(p)c R (k)f 



where 
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and TRi(m) = T R ^(kF,m)TLi(kF,m) Based on 

Ht, the tunneling current operator J is given as follows: 



J = eiY, [c f R (k)c L (p)fl - c[(p)c R (k)f 

k,p 



(8) 



We now compute the expectation value of current op- 
erator via a conventional first-order S-matrix expansion: 



J(t) 



dt'([J(t),H T (t')}). 



(9) 



The new aspect of our tunneling theory is the vertex 
operator T which contains all of many-body effects of 
the exciton condensate. Eq.@ can be evaluated further 
using the Feynman diagram depicted in Fig|5J 

/°° dp 
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x [f(s)-f(e + eV)} 
= 4Tre 2 D R D L \(f)\ 2 V 



FIG. 3: Normalized expectation value of condensate order 
parameter (S x ) 2 ■ Dashed lines indicate the upper and lower 
bound for an estimate of the thermodynamic limit of (S x ) 2 
as a function of cL/Ib- N is the total number of electrons. 



where A R (At) is the spectral function of the right (left) 
lead, /(e) is the usual Fermi-Dirac distribution function, 
and D R (Dl) is the density of states at the Fermi sur- 
face of right (left) lead, ft is clear from Ea. Hl(J|) that 
there is no DC Josephson effect because the conductance 
G (= dl/dV oc |(T)| 2 ) is finite. However, the interlayer 
tunneling current is zero unless there is a phase coher- 
ence: (T) ^ 0. Remember that (T) measures a phase 
coherence between states with different values of inter- 
layer number difference, AT re i, because T cx c^Cmi and 
therefore changes 7V rc i by two. So, unless the ground 
state is a coherent linear combination of states with var- 
ious iVrei, (T) is zero, and so is the tunneling current. As 
mentioned before, this is similar to the phase coherence 
between different number eigenstates in superconductiv- 
ity, which is responsible for the Josephson effect. In this 
sense, interlayer tunneling conductance is related to the 
Josephson effect. However, we emphasize that the con- 
ductance should be finite even at zero temperature and 
there is no direct analogy with the Josephson effect. We 
now compute the interlayer tunneling conductance as a 
function of d/l R - In particular, we will be interested in 
normalized conductance since the absolute scale of con- 
ductance is sensitive to sample-specific details such as 
D Rl D L and T RL . 

In essence, we compute |(T}| 2 which can be further 
reduced as follows: 



is*) 2 



(11) 



(10) 



where N is the total number of electrons, and we have 
used the fact that (c^Cmi) is independent of m and is 

equal to (S x )/N. S x [= E m ( c LT c »U + c Li c ™t)/ 2 ] is tlic 
order parameter of exciton condensation, and it can also 
be interpreted as the pseudospin magnetization in the 
x direction. Since ^2 m Tsj J (m)/N does not depend on 
d/l R , the interlayer distance dependence of conductance 
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FIG. 4: Coulomb self-energy of a quasiparticle. 




FIG. 5: Normalized interlayer tunneling conductance peak 
as a function of interlayer distance in comparison with ex- 
perimental data from Ref. lj. We define the normalized con- 
ductance as conductance divided by its maximum value as a 
function of d/ls- Two theoretical curves are obtained from 
the upper and lower bound of thermodynamic estimate in 
FigEl 



is solely determined by {S x ) 2 / A 2 ^. 

In FigO we plot (S x ) 2 as a function of d/ls which 
is computed via exact diagonalization of finite systems 
with various particle numbers in torus geometry. When 
computing (S x ) in finite systems, it is very important to 
take into account fundamental fluctuations in iV re i; the 
true ground state is a coherent, linear combination of 
states with various iV re i 0, 0] . Though estimating the 
accurate thermodynamic limit of (S x ) 2 is difficult, it is 
reasonable to argue that the true thermodynamic limit 
lies between two dashed lines in Fig|3J 

FigEl plots the Coulomb self-energy of a quasiparticle, 
Ac /(e 2 /els), as a function of d/ls which is determined 
in exact diagonalization studies by computing the energy 
gap of particle-hole-pair excitation with the largest mo- 
mentum and taking half of its value. For comparison, we 
also plot the self-energy in the Hartree-Fock approxima- 
tion |l8| which tends to overestimate Ac- 

Finally, in FigO we compare our estimate of normal- 
ized interlayer tunneling conductance near zero bias, i.e. 



{S x ) 2 /A^, with experimental data of Spielman et al. Q]. 
We define the normalized conductance as conductance 
divided by its maximum value as a function of d/ls- 
Two dashed lines in Fig[S] correspond to the upper and 
lower bound of estimated thermodynamic limits of (S x ) 2 
in Fig|3J Considering simplifications used in our the- 
ory such as omission of finite thickness effect, we find 
our theory to be in good agreement with experiments. 
In addition to further comparison with experiments in 
the regime d/ls > 1-2, it will be very interesting to see 
whether our prediction of decrease of conductance peak 
for d/ls < 1-2 is consistent with future experiments. Re- 
member that decrease in conductance peak at small d/ls 
is due to increase in energy gap to put electrons into bi- 
layer systems while the pseudospin magnetization is sat- 
urated. We would like to emphasize that, once normal- 
ized, our theoretical estimate of conductance peak does 
not have any fitting parameter. 

We have shown by means of Eq. (|10|l and that in 
exciton condensate the interlayer tunneling conductance 
at small bias is finite, but strongly enhanced. However, 
we did not show why the conductance should be sharply 
peaked near zero bias, which we will explain now. Once 
the interlayer current is driven by an external electromo- 
tive force, it should physically flow through the bilayer 
system since otherwise there is no steady state. Exciton 
condensates accomplish this by adjusting their interlayer 
phase difference <fi to sustain the externally driven cur- 
rent, which is again easy to understand in terms of the 
ground state wavefunction at d/ls — > 0: 



IV>mW> = II(4 



s*4i)|0>, 



(12) 



which carries a net internal current within bilayer sys- 
tem equal to eAsAS y sm ^ E3- Then, there should be 
a critical current at <j) = tt/2 which is the maximum cur- 
rent allowed without breaking phase coherence. There- 
fore, for sufficiently large voltage bias, coherent inter- 
layer currents should be cut off and become constant 
as a function of bias voltage, once they reach the crit- 
ical value controlled by single-particle interlayer tunnel- 
ing gap AgAs- The conductance associated with coher- 
ent tunneling, therefore, should be zero after the criti- 
cal voltage and is strongly enhanced only near zero bias. 
Consequently, the width of conductance peak is propor- 
tional to very small, but finite AgASi while the propor- 
tionality constant strongly depends sample-specific de- 
tails such as the density of states of leads. It is, however, 
encouraging to find that typical width of conductance 
peak (~ 10 — 100[ieV) is roughly in the same order as 
Asas Gill ■ The above argument is valid for general d/ls 
when there is phase coherence. 

Until now, we have studied the interlayer tunneling 
conductance in a single bilayer system, which, we showed, 
is not the exact analog of Josephson effect. We now con- 
clude by proposing a much more direct analog with the 
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Josephson effect. Consider a pair of bilayer systems, say 
A and B (four layers altogether), separated by a lat- 
eral tunneling barrier. Then, put an interlayer current 
through the top and bottom layer of, say, bilayer system 
A, in which way a non-zero interlayer phase difference 
is induced in bilayer system A while the system B has 
none. We predict then that there will be two counter- 
flowing currents: one between two, top layers of system 
A and B, and the other between bottom layers. The net 
current will be zero, but it may be possible to measure 
these two currents individually. 
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